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Abstract. Flux qubits, small superconducting loops interrupted by Josephson 
junctions, are successful realizations of quantum coherence for macroscopic variables. 
Superconductivity in these loops is carried by ~ 10^ - 10^° electrons, which has 
been interpreted as suggesting that coherent superpositions of such current states 
are macroscopic superpositions analogous to the state of Schrodinger's cat. We 
provide a full microscopic analysis of such qubits, from which the macroscopic 
quantum description can be derived. This reveals that the number of microscopic 
constituents participating in superposition states for experimentally accessible flux 
qubits is surprisingly but not trivially small. The combination of this relatively small 
size with large differences between macroscopic observables in the two branches is seen 
to result from the Fermi statistics of the electrons and the large disparity between the 
values of superfluid and Fermi velocity in these systems. 
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1. Introduction 

Schrodinger's cat paradox [T] emphasizes how basic quantum concepts such as 
superposition that are routinely apphed and accepted in the microscopic description of 
matter appear to contradict basic human experience when augmented to a macroscopic 
scale. Quantum mechanics does not itself provide any intrinsic size limitation, i.e., it 
does not predict a critical number of particles at which such superpositions would be 
impossible. Yet there still remain questions as to whether there are emergent limitations 
at some length scale between the micro- and macroscopic, or even whether extraneous 
factors limit the size of superpositions. This question is fundamentally unresolved [2] . 
Hints towards an answer would be given by experiments that realize truly macroscopic 
superpositions, which would provide evidence against macroscopic realism [3]. There 
are numerous experiments attempting to produce such macroscopic superposition states 
in a variety of different physical systems ^ El El [71 El [9] . Yet determining the actual 
size of these states (sometimes referred to as Schrodinger's 'cattiness' [2]) is a non- 
trivial theoretical question that needs to be solved in tandem with the experimental 
realizations. This includes three aspects. First, one must ascertain or specify which 
types of degrees of freedom participate in the superposition, and which are irrelevant. 
For example, trapped ions in a superposition state of internal degrees of freedom interact 
with the motional degrees of freedom in the ion trap, yet these would not be included in 
a superposition size measure for the internal states. More generally, degrees of freedom 
that are not accessible in a particular experimental realization, e.g., because of too 
high energy, are not included. Thus in a low-temperature experiment one would not 
determine the superposition size in terms of quarks. Second, once it is agreed which 
elementary constituents are involved, one needs to determine what is the actual number 
of particles or modes that participate in the superposition, in the sense of being in 
distinguishably different states in the branches of the superposition. The latter can be a 
highly non-trivial calculation for an interacting quantum system. Third, it is important 
to then assess this number in the context of the observable and controllable physical 
parameters of the superposition state. 

One attractive candidate system in this program is the superconducting flux 
qubit. Flux qubits are composed of superconducting loops of between one and 
many micrometres, that contain one or more Josephson junctions. They can realize 
superpositions of states of a macroscopic electrodynamic variable, the circulating current 
and its concomitant mag netic flux, attributable to 10^-10^° electrons |2l [3 El [ID]. 
Superpositions of the form |0) + |0) [Zl [H] and coherent oscillations between these 
states have both been demonstrated [HI |12]. However, even though many particles 
are involved and the fluxes are macroscopically distinct in the two branches of the 
superposition, it does not necessarily follow that the actual size of the state, i.e., the 
number of particles that are in superposition in the sense of being in different states in 
the two branches, is macroscopic as well. Furthermore, this size can not be determined 
from the engineering of these experiments, for which the macroscopic electrodynamic 
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variables are sufficient [3, El HOI IHl [E] • Finally, a microscopic calculation has to carefully 
take into account that the superconducting states is composed of (paired) electrons. 
Because of their Fermion nature, these always occupy a finite volume in momentum 
space (unlike Bosons). Moreover, they are indistinguishable. It is thus important to 
devise a measure that does not assume distinguishable electrons and that correctly 
takes both the indistinguishability and Fermion statistics of the electrons into account. 
Theoretical estimates have come to vastly different proposals for this size, reflecting 
their different underlying assumptions: i) counting all electrons of Cooper pairs in the 
supercurrent gives a size of 10^-10^° [21 [3 El [13], which is at least mesoscopic, if not 
approaching macroscopic sizes, ii) analyzing the number of electrodynamic charge states 
within the macroscopic, circuit Hamiltonian approach leads to a trivial effective size of 
1-2 [H]. However, neither of these approaches is fullly microscopic in the sense of 
characterizing the interacting electron system in the superconducting loop. 

Here we place a bound on a broad range of microscopic measures of superposition 
size by asking the following direct question: how many electrons behave distinguishably 
differently in the two branches |0) and |0)? We first outline the microscopic analysis 
of the superposition state that is necessary to estimate this quantity and then we 
present the bound for electrons in superpositions of flux states. This bound derives 
from an operational measure for superposition size that is determined by the degree 
of distinguishability of the two branches by an n-particle measurement [15]. We will 
show that for all experimentally realistic flux qubits, this bound results in an estimate 
of superposition state size that is considerably smaller than the number of electrons 
carrying the supercurrent but that is also significantly larger than the trivial size of 
1-2 estimated from the macroscopic description. Our microscopic analysis reveals that 
the Fermi statistics play a critical role in reducing the number of electronic degrees of 
freedom that participate in the quantum superposition to well below simple estimates 
of the number of current carrying electrons. We further discuss the questions raised 
by such an estimate of size in terms of microscopic constituents, for a system whose 
macroscopic behavior can also be fully described by a single collective variable. 

2. Microscopic analysis 

In this Section we outline the analysis performed to connect the full microscopic 
description of the system, starting from the Hamiltonian of Bardeen, Cooper and 
Schrieffer (BCS), to its macroscopic state and show how to compute microscopic 
quantities such as density matrices and correlation functions in that state. We consider 
a flux loop containing a single tunnel junction and employ the functional integral 
formalism for superconducting tunnel junctions that was ffist developed in Ref. [16] . 
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2.1. Hamiltonian 

The BCS Hamiltonian for the electrons in the bulk of the superconductor has the 
following form (in Gaussian cgs units): 



BCS 



2m 



V 
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where ipa is the electron field operator, e is the (negative) electron charge, and the sum 
is over spin indices a = ±1/2. is the electromagnetic scalar potential The quartic 
term is the effective Cooper pairing electron-electron interaction, and includes both 
electric forces and phonon interactions||] The free electromagnetic field is governed by 
the following Hamiltonian: 



EM 



d=^r 
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] = y^^r^ (^-^Aind - V<^ind^ +(VxAind) 



(2) 



where h is the magnetic field and A^d = A — Acxt- Only the fields Ejnd = E — Ecxt 
and hind = h — hcxt induced by the electron dynamics are connected with the electron 
Hamiltonian, and therefore we subtract the external fields from the total field in Hem- 
Note that the integral in Eq. ^ is over all space. 

An effective description can now be developed as follows: 

(i) We assume that there are no free electric fields in the bulk of the superconductor, 
so that the only terms in this region will be H^cs and the term (V x Aind)^ of 

-f^EM- 

(ii) Both close to and inside the Josephson junction, any induced electric fields across 
the junction and their interaction with the surface electrons can be described by 
an effective capacitive term ^CV"^, where V is the line integral of — E = ^A + Vip 
across the junction. 

(iii) Outside the conductor, the energy of the magnetic field can be re-expressed as an 
inductive term, — $cxt)^) where $ is the flux enclosed by the superconducting 
loop. 

(iv) The dynamics of electrons close to the junction due to tunneling across the junction 
can be expressed as an effective tunneling Hamiltonian, 



Ht = J2jd'rj dV^t(r)Trr'^.(r') 



(3) 



where Trr' is a tunneling amplitude which is non-zero only when r and r' are close 
to and on opposite sides of the junction. 



I This implies that the scalar potential ip only includes potentials due to external fields, not generated 
by electron dynamics in the bulk of the superconductor. 
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The effective Hamiftonian of tlie total system, bulk superconductor and tunneling 
junction, is then 

H = H^cs + i^T + ^ / (V X Aind)' + (-^ - '^'cxt)' + (4) 

OTT J ZL L 

where i^T is only relevant close to a junction and the integral over the vector potential 
is now taken only over the interior of the superconductor. 



2.2. Expectation Values 

The electron field operators can be written in terms of creation and annihilation 
operators and Ck,o- for plane wave modes with wave vector k as 

Mr) = J d^kc^,^e-'''-\ (5) 

Our goal will be to calculate arbitrary correlation functions in the grand-canonical 
ensemble between electron field operators. This will then also allow us to calculate 
expectation values of any operator O that is a function of electron creation or 
annihilation operators. We shall be specifically interested in mode occupation numbers, 
both for single electron modes and for Cooper pair modes. The former are given by 
2-point 1 mode correlation functions of the form ^cj^o-^k.o-^ and the latter by 4-point 2 

mode correlation functions of the form (^c]^ ^c^_y. ^C-k^iCk^^"^ . To achieve these goals, we 
use the finite-temperature functional integral formalism with imaginary time-variable r 
and electron variables anti-periodic in r, having period P = 1/kT [T7] All expectation 
values calculated here will be restricted to the pure ground state by taking the limit 
/3 -> oo. 

The desired expectation values of an operator O are given by the functional integral 
expression [17]: 



with action functional S[^}j\ given bj|§] 



d 



+ H 



where /i is the chemical potential and H is given by Eq. Q. We define the generating 
functional Z [^] [IZl [HI [E] : 

§ Going to imaginary time means that an extra factor of i is picked up by if, which formally is the 
time-component of the four- vector , so that the potential term of i?BCS becomes —ieip'^ip'tp^. 
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One can now eliminate the quartic term in Hbcs by introducing an extra integral 
over an auxiliary field and applying the Hubbard-Stratonovich transformation [20 
fT7] . This leads to S[ip, A, cp] in Eq. ([t]) being replaced by the effective action: 



+ ^ ir.r^M'^ + Ae-*<^ ij^^ij,) +^A^ + Ht + Hem } 

where the auxiliary field Ae**^, with A real-valued, has been introduced, and where the 
functional integral in Eq. ([s]) now also runs over all configurations of the fields A and 



, (9) 



For convenience, it is usual to make a transformation ipo- — ?• ipa^'^'^^'^ to remove the 
factors of e^*''^ from Eq. ([9|. The effect of this is that the operator V is replaced by 
V — iV(f)/2, which can then be combined with the vector potential A to define the 
gauge- invariant superfluid velocity vs: 



^ /^^ , 2e , 
Vs = V0 + — A 

Im V no 



(10) 



Here m the single electron mass, e the (signed) single electron charge, and pe the 
superfluid electron density, i.e., the effective superconducting electron density [22]. vs 
is equal to the mean velocity of superconducting electrons in the system and is related 
to the superconducting electron current density via j = epe'Vs- Note that the density 
of superconducting electrons, pe, is twice the density of Cooper pairs, ps (commonly 
referrred to as the pair superfluid density), while vg is equivalent to the center of mass 
velocity of Cooper pairs. 

Further introducing the Nambu space [23] in which the electron field is expressed 
as the vector 



transforms the electronic part of the effective action into a proper quadratic form: 
Sei = J drj dV*t(^,r)G-i(r,r;r',r')*(r,r'), 



with 



where 



^-1 

•^bulk 



-1 



G 



h? im 
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(12) 
(13) 



p — ICif 



Act, 



and the tunneling matrix T is given by 



Trr' 



Trr' 



5{t-t'). 



(5(r - r')(5(r - r') (14) 



(15) 
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The matrix G^^^y. is seen to be the inverse of the Green's function for the field equations 
for the electron field in the bulk of the superconductor. 

The generating functional Eq. dsj) can be simplified by integrating over the electron 



field since the electronic part of the action Sei, Eq. (12), is now quadratic in ip. This 
results in a Gaussian integral which yields 



The functional Eq. (16) gets its largest contributions from those field configurations for 
which S[A,ip, A,(j)\ has a maximum, i.e., where the functional derivatives of S with 
respect to each of the fields is zero, i.e. a saddle point. For fields where the second 
derivative is very large, S will get its only significant contribution from the saddle point 
values, i.e., the fields will behave classically. As we described below, this turns out to 
be the case for the fields A, ip and A but the second derivative of the effective action 
with respect to the phase is not necessarily macroscopic. Hence we cannot (at the 
moment) fix to a single classical value, and this remains a quantum variable over 
which the functional integral is to be carried out. 

We now summarize the results of the saddle point analysis for A, ip, and A. For A 
we find 

X V X (A - A„) = j = (V* + ^a) . ep.v. (17) 

The identification of the far left-hand side with the current density j is done using 
Maxwell's equations (Ampere's law). For ip we find 

= 2eip, (18) 

which we recognize as the AC Josephson equation. For A, since in the situation we 
are interested in, the magnetic and electric fields will be much weaker than any critical 
values or energy scales and slowly varying, we may expand trlogG"^ perturbatively. 
This leads to the following saddle point equation: 



ppg 



/3 



dE , tanh - + = 1, (19) 



where pp is the density of states at the Fermi surface, ujb is the Debye energy, and 
E denotes the difference in energy between a given energy level and the Fermi energy. 
Note that this is equivalent to the self-consistency equation for the energy gap A in 
BCS theory. As all currents, fields, and temperatures in the system are far below their 
critical values, we can thus replace A with its standard BCS value. 

We are left with as the only non-classical variable. Expanding tr log close to 
the junction perturbatively in Trr' leads to the further simplification 

Z\i\ = /p0e-^«''^[^^*'-^^™^^+^(*-*-*)'+^''^^'^^], (20) 
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where the flux $ is related to the phase variable through the flux regular flux 
quantization relation 

^=i-^o, (21) 

and the Josephson energy Ej is an effective quantity whose form depends on the details 
of the tunneling matrix elements (a simplified form is given in Ref. [21] )• 

We now obtain correlation functions (expectation values of products of creation 
and annihilation operators) by taking time-ordered functional derivatives of Z [C,] with 
respect to This yields products of Green's functions which are linear combinations 
of products of the equal-time Gorkov Green's functions US]- We note that, as 
derived also in Ref. [16], the zero temperature limit of these Green's function products 
yields dynamics for the order parameter that is equivalent to that of a particle in an 
effective potential Ej cos ~ ~ 'I'cxt)^, which is precisely the description used in 
the macroscopic theory of flux qubits [12j. 

The full Green's function, G, is written as a sum of the zeroth order bulk Green's 
function and perturbative contributions from the superfluid flow and junction tunneling: 

G = (Gb,\k + Trr'y'= Go + 5Gvs + SGt. (22) 

In order to compare states with different current distributions, the Green's functions 
must be expressed in terms of velocity relative to a stationary laboratory frame, i.e., 
using modes that are eigenfunctions of the laboratory frame momentum operator 

/iq = —ihV + mvs, (23) 

rather than of /ik = —ihV. These bases can be viewed as related by a Doppler shift. 
Given that |vs| is smaller than the critical velocity, the energetics can be assumed not 
to change between frames [25]. As we describe in Section |4| the correction 5G^^ is 
evaluated by perturbatively expanding Go(q, q') in terms of the superfluid velocity vs. 



This yields [IHl [19] (see also discussion in Section 4.1): 

^^""^ o^c.. , A.N 3/2 



5G(q, q') = -(27r)35(q - q') A (Ao", - J . " (24) 



2(E2 + A2)- 

where (t^^cTz are Pauli spin matrices in the Nambu space. The first order contribution 
to the Green's function from the tunneling matrix, 5Gt, is obtained as [HI [19] (see also 



discussion in Section 4.2): 



jim SGA^L = q^) = ^ ( (\ + ,-.A.)^ ^ (\ _ ,.A.| 



(25) 



Here Qk = V -^l ~^ where Ek is the electron kinetic energy, Ek = — Hf, with /xp 
the Fermi energy, and A7 = A0 — ^ J^^ A • dr is the gauge invariant phase difference 
across the tunnel junction. 
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3. Superposition size measures 

The size of a superposition state of a many-particle system is not a unique physical 
observable defined by its experimental measurement procedure, leaving many ways how 
to define superposition size. A naive way to define such a size would be to compare the 
absolute value of the difference in some suitable physical observable between the two 
branches to some characteristic atomic scale for that observable. However, as pointed 
out in Refs. |2l EZj, this approach of estimating an "extensive difference" between the 
branches is too simplistic. For example, the superposition of a single neutron going 
two ways through an interferometer involves a huge difference in both center-of-mass 
position and angular momentum around the center of the interferometer compared to 
atomic scales, yet still involves only a single neutron and would clearly not be considered 
a macroscopic superposition state [2]. Similarly, for superposition states of flow, e.g., 
in flux qubits, merely comparing the difference in magnetic moment A/i between the 
branches |0) and |0) to the magnitude of the Bohr magneton fiB can yield the same 
value A/i = n^s for very different n particle states. For example, neglecting electron 
indistinguishability, a state in which n electrons are in orthogonal states in the two 
branches, each having the same angular momentum difference AL = h, results in the 
same A/x as a state in which a single high energy electron carries the equivalent total 
angular momentum difference AL = nh while all other electrons behave identically 
in both branches. Only the former case involves n electrons behaving distinguishably 
differently in the two branches and would be acceptable as a superposition state of 
size n: in the latter case only a single electron behaves differently and this is clearly a 
superposition state of size unity. However, for electronic states of flow, we also have to 
take into account the fact that electrons are themselves indistinguishable, with Fermi 
statistics. This further complicates the analysis, requiring an accounting of how many 
electronic modes behave differently in the two branches, i.e., have different occupation 
numbers, rather than how many individual particles behave differently. This example 
suggests that in addition to differences between macroscopic variables, one needs to 
i) consider also the number of microscopic constituents that are behaving differently 
in the two branches when defining the effective size of a many-body superposition 
state, and ii) take the fundamental particle statistics into account when defining these 
microscopic constituents. Since electrons satisfy Fermi statistics, the correct microscopic 
constituents for analysis of superposition states in flux qubits are modes of electron 
momentum, occupation of which can be different in the two branches. 

Several possible definitions have been put forward for a superposition size measure 
based on analysis of microscopic degrees of freedom [HI ETJ [15] . While all of these are 
hard to calculate for complex many-body states such as the electron states constituting 
a flux qubit, the distinguishability measure developed by two of us in Ref. [15] admits an 
upper bound that can be calculated using the microscopic functional integral formalism 
described above. The distinguishability measure is an operational measure that asks 
what is the largest number of subsets of elementary constitutents such that measuring 
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all constituents in any one subset causes the superposition to collapse to one branch 
with some specified probability. This problem is equivalent to that of determining the 
minimum number of microscopic constituents that have to be measured in order to 
distinguish the two branches to some specified precision. Ref. [T5| showed that this 
leads to a superposition size of N/rimin, where n^m is the first value for which the 
probability of successfully distinguishing the branches is larger than 1 — 5, where 6 is 
the desired precision. We adapt this here to consideration of the number of electron 
modes n that need to be measured in the mode electronic system. For simplicity we 
give the analysis here in terms of single electron modes. It is important to note that 
the single electron modes combine to make Cooper pairs rather than excitations. As we 
show in more detail below, identical results are obtained if Cooper pair modes are used 
in our analysis. However, we prefer to make the primary analysis with single electron 
modes since these highlight the underlying Fermi statistics, which as noted earlier are 
also refiected in the fact that the BCS state occupies a finite volume in momentum 
space. For small differences in mode occupations 



^(n) ^(n) 

Pa -Pb 



1=1 







l-iV„ 



/ 1 -AN„ 



(26) 



Averaging the probability of inferring the correct branch by measuring n selected modes 



P. 



1 1 



Pa 



Pb 



1 



(27) 



over all possible choices of n out of JV modes leads to 



P. 



tot 



(28) 



1 n 

where the second term on the right hand side is an upper bound, derived without 
regard to multiple mode occupancy. From this a lower bound on the limiting value 
of n for success is extracted as n = A/'/AA^tot, resulting in an upper bound for the 
superposition size of AA^tot, i-e., the total difference in occupation number between the 
two branches. This estimate places an upper bound on any such microscopic measure of 
superposition size by asking how many electrons on average are in a different mode in 
one branch relative to the other branch. Any measure that gives a larger number than 
this must be counting electrons that are in identical states in the two branches and are 
therefore not actively contributing to the superposition. Since our measure is based on a 
measurement, it is desirable (although, just as in a gedanken experiment, not essential 
for theoretical evaluation of the superposition size) that the assumed experiment is, 
in principle, feasible. Thus it is essential that the indistinguishability of electrons be 
incorporated in the superposition size measure. ANtoti the total change in occupation 
numbers of all electron modes in the system is indeed the only meaningful indicator of 
how many electrons are affected when passing from one branch of the superposition to 
the other, when the indistinguishability of electrons is taken into account. We return 
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to this in Section [5] when discussing the role of Fermi statistics in determining the 



numerically obtained values for superposition states of flux qubits (Section 4.1). 

For our calculation based on the BCS theory |28], as outlined above, the natural 
choice of single electron modes are momentum eigenstates. We now show that this 
choice of single-particle basis also maximizes the occupation number differences between 
the two branches, 6nq = (0| c^^^Cq^o- 10) — (0| c]j o.Cq,o- 10), which ensures that AA^tot 
is a mode-independent and hence true upper bound on the superposition size. Sn^ 
can be viewed as the diagonal elements of a matrix D with elements -Dqcqv = 
(0| c^^Cq/,^' |0) — (0| c]j_^Cq/^^/ |0). AA^tot = Z^k l^'^kl IS then the sum of the absolute 
values of the diagonal elements of D, i.e., AA^tot = tr ||D||. It can readily be shown that 
this number is maximal in the basis in which D is diagonal [181 HH] • For a superconductor 
we have (cl^^c^iy'^ oc 6(k — k')6aa', so D is diagonal in a momentum basis. Therefore 
the occupation number differences must be evaluated in a momentum basis and not in a 
position basis, in order to correctly evaluate our upper bound to the superposition size. 

Since we use BCS theory, the effect of Cooper pairing on the electron dynamics is 
automatically included, but to confirm this explicitly we shall also define two-electron 
Cooper pair modes and show that using these modes yields the same result for branch 
occupation number difference as that obtained using the single electron modes. 

4. The number of electrons changing modes between branches 
4.I. Loop current contribution 

Superposition states in flux qubits can be described as superpositions of clockwise 
and counterclockwise net circulating currents of the form (|0) + |0)). While 
these currents have intrinsic quantum fluctuations, these fluctuations do not affect 
the quantitative estimates made here which make use of the experimentally measured 
average currents (and hence implicitly an averaged superfluid phase parameter, since the 
latter is linearly related to the current). The saddlepoint solution of our path integral 



equations described in Section 2.2 is then characterized by the following assumptions. 
For qubits with thickness much smaller than the penetration depth A, we can take the 
flow to be approximately uniform in the lateral dimension, i.e., the thickness in a quasi- 
planar geometry. The current flow is conveniently characterized by the complex phase (p 
of the superconducting order parameter A, and on the electromagnetic vector potential 
A, combined in the gauge-invariant combination that defines the superfluid velocity, 



Eq. (10). From the structure of the collective variables at the saddlepoint, the Green's 
functions and hence the occupation numbers are obtained by Fourier transforming and 
taking imaginary time-ordered functional derivatives of the generating functional Z [^] . 
We now summarize the main points of this calculation. Full details will be given 
elsewhere [19] . 

Recognizing that as long as both the externally applied and the internally generated 
magnetic fields are weak and the cooling is sufficiently adiabatic to avoid vortex 
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formation, the entire electron system may be decomposed in terms of the single electron 
modes. The 2-point, 1 mode correlation function is obtained as 



{^\ci^c^,\^) = hm (V^:(k,r)V^.(k,0)) (29a) 

= lim . Z ^] 29b 

= - lim hG^Jk,T';k,T), (29c) 

where G„fj is a matrix element of G. As noted earlier, in order to compare properties 
of states with different current flow, in particular the two branches of a flux qubit 
superposition with left and right circulating currents, the correlation functions must be 
expressed in terms of the laboratory frame wave vector q = —iV + mvs/^. The loop 
current contribution to the single mode correlation functions is then obtained as 

(^Ict^Cq^l^) = -/i^hm a,(q-^,0;q-^,r). (30) 

Provided that the current in the superconductor is small compared to the critical current 
-^c,buik, a condition which is satisfied in the flux qubit experiments [HIIZIIH]! the superfluid 
velocity vs will be a small perturbative quantity and calculations can be carried out to 
first order in |vg|/fcrit = |vs|fF"^/^5 where vp is the Fermi velocity. Expanding Go-o- 
leads then to 

clc^^ = ^ - + ■ ' (31) 

where (vs) is the mean superfluid velocity averaged over the quantum state of the 
system. The average occupation number of a single electron mode is given by nq = 
(^q,o-^q,o-) (^q is independent of a for any reahstic magnetic field strength). Thus 
Eq. (31) provides the contribution from Gq and SG^s to the mode occupation number. 



The difference in occupation number 5nq of a mode (q, a) between the two circulating 
current superposition branches, to first order in vs, is then obtained as 

^'^q = (01 c^^^cq,. 10) - (01 c^,,cq,. |0) = ^ " '^(vs) • (32) 

Before evaluating this difference for the recent flux qubit experiments, we 
summarize the corresponding analysis in terms of Cooper pair modes. A Cooper pair 
mode accommodating two Cooper paired electrons with intrinsic momenta and spins 
(k, 'I') and (— k, J,) is defined by the two-electron creation operator Cl. = cj^i-ct^kj,- '^^^ 
Cooper pair mode correlation function (cl-C^^ is equal to the two-point correlation 

function for two single electron modes, (^cj^^c^,^ j^c_k,4,Ck,t^ and the branch occupation 
number difference for Cooper pair modes is equal to 

5^fk,-k = (0| 4tC-HC-HCkt |0) - (0| 4t^-HC-k4Ckt |0) . (33) 



Electronic structure of superposition states in flux qubits. 



13 



Analysis of the corresponding correlation functions in the laboratory frame is 
complicated by the fact that two laboratory frame velocity-modes q = k + mw^/h 
and q' = k' + m\^/h with opposite spins and with wave vectors k and k' relative 
to vs will be coupled by the Cooper pair coupling only if k' = — k. Hence q' = 
— q + 2mvs//i and so whether two modes are coupled or not will in principle depend 
on the superfluid velocity. This can lead to a change in mode definitions between 
the two circulating current states that complicates the calculation of mode occupation 
number difference. However, since for the experimental systems nw^/h is much smaller 
than the momentum-space difference between distinct modes, this mismatch can be 
treated perturbatively. Detailed analysis shows that if the laboratory frame modes 
are not perfectly correlated, i.e., q' 7^ — q + 2m\a,/h for the value of vs in either of 
the branches, there can be only a second-order dependence on vs, while if they are 
perfectly correlated, i.e., q' = — q + 2mvs//i in one branch, a first order contribution 
can result [HI [19]. Furthermore, the maximum value of this first order contribution is 
exactly the same as that obtained from the single electron mode occupation number 
difference, i.e., (5A/'q,q' = 5nq5(q + q' — 2mvs/^). Thus, as might have been expected, 
the two participating single electron modes are perfectly correlated. More importantly 
for the present discussion, this result imphes that analysis of the branch occupation 
number difference of Cooper pair (two-electron) modes yields exactly the same effective 
superposition size as that obtained from analysis of the branch occupation number 
difference of single electron modes. 

We now use the single electron mode analysis to derive an expression for the branch 
occupation number difference and hence an upper bound on the effective superposition 
size in terms of experimentally measured quantities. We start by finding the total 
number of electrons changing mode in a local spatial region, which is obtained by 



summing Eq. (32) over all modes q: 

An{r) = 2vr y dqq'^pq^^^ J d(cos6') |5vs(r)| cos 6' 



9 

TT/Ci |(5vs(r)|. 



(34) 



Here 6 is the angle between q and 5vs, Pq is the (unknown) density of modes per unit 
volume and we have denoted the integral of A'^hq^ /2Q^ over all modes q by /Ci. In 
addition, to avoid double-counting electrons, i.e., counting them both as they leave one 
mode and enter another, we only sum over the modes for which 5nq > 0. Using the 
relation between current and velocity, we may derive a related expression for 5j(r), the 
local current density: 

5j(r) = e [ d'qp,\5n{ci)\^. (35) 



m 



From Eq. (32) it is evident that components of q perpendicular to (5vs give zero 
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contribution, so we can replace q by its parallel component q cos 9 to arrive at 

5j(r) = 27re / dqq^pq—— |(5vs(r)| / d(cos6') cos^ 6* 

= ye/Ca |5vs(r)| , 

where /C2 denotes the integral of A^/i^g^/2mf2^ over all modes q. To evaluate the 
integrals JCi and IC2, we note that the denominator in Eq. (32) strongly suppresses 
modes away from the Fermi surface, i.e., Sn^ is non-negligible only for modes close to 
this, as expected for modes participating in Cooper pairing, pg and q may then be 
replaced by their values at the Fermi surface, pp and qp = y/2mp/h, respectively, while 
the q dependence of Qq must be maintained since this varies significantly over the range 
of q. The integrals then yield /Ci = pY2m?p/h?' and /C2 = {hq^ /m)lCi, which allows the 
unknown pq to be eliminated, to obtain the local occupation number difference in terms 
of the local current difference between the two branches: 

Mr) = ipM. (37, 

Finally, integrating over the entire volume of the superconductor yields the total branch 
occupation number difference AA^tot, the total number of electrons that are in different 
modes in the two branches, as 



3/ 

AiVtot = ^ (3J 



where L is the total length of the main superconducting loop of the flux qubit and 
5/p is the experimentally measured difference in "persistent current" [7j between the 
superposed branch states. We note that this derivation assumes that the current 
distribution is homogenous on scales less than or equal to the Fermi wavelength but 
not necessarily beyond this length scale. Also, while the use of experimental persistent 
current values implicitly incorporates an average over quantum fluctuations in the 
superconducting phase parameter 0, explicit measurement of current fluctuations could 
allow a higher order analysis. 

We have evaluated the occupation number difference bound on the effective 
superposition size for all reported experimental demonstrations of flux state 
superpositions to date. The relevant experimental parameters and corresponding values 
of AA^tot, the effective superposition size, are listed in Table [T] for the three recent 
experiments [El [TJ |9]. The largest numbers are found for the SUNY experiment [8] 
carried out between excited states of a single- junction RF-SQUID configuration, while 
the Delft p'J and Berkeley experiments were both made with three-junction flux 
qubits that generated a superposition of degenerate ground states. The latter are very 
different in their geometric size. We also list in Table [T] the corresponding values for the 
difference in the two macroscopic observables current and magnetic moment between 
the two branches, respectively 5Ip and 5p = ASIp, where A is the area enclosed by the 
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Experiment 


Material 


Vf 


L 


5/p 


Afj, 




SUNY 


Nb 


1.37 X 10*^ m/s 


560 fim 


2-3 /iA 


5.5-8.3 X IOVb 


3800-5750 


Delft 


Al 


2.02 X 10^ m/s 


20 /im 


900 nA 


2.4 X IOVb 


42 


Berkeley 


Al 


2.02 X 10^ m/s 


183 /im 


292 nA 


4.23 X IOVb 


124 



Table 1. Parameters and effective superposition sizes for current superposition states 
produced at SUNY [Sj, Delft [7] and Berkeley wp is the Fermi velocity, L the 
length of the superconducting loop, 6Ip the measured difference in persistent current 
between the two branches and (5/i — ASIp is the different in magnetic moment, where 
A is the area enclosed by the loop. AA^tot is the effective superposition size, i.e., the 
total number of electrons participating in the superposition state. 



superconducting loop. Before discussing these results, we first briefly summarize the 
possible contributions from tunneling through the junction and the effect of scattering 
from impurities. 

4-2. Tunneling contribution 

The estimates in Table [T] include only the effect of different circulating currents in the 
two branches of the superposition. This is the only contribution in the bulk of the 
superconductor, far from any junctions and where the tunneling contribution to the 
Green's function G is negligible. However close to the junction, 6Gt is not negligible 
and electron tunneling through the junction can then also contribute to the difference in 
mode occupation in the two branches of the superposition. Physically, this corresponds 
to the fact that close to the junction up to a depth of order = vp/A, the Cooper 
pairs will be in superposition between both sides of the junction, rather than being 
described by rigid, directional superflow. This gives rise to an additional contribution 
to the superposition size when the modes are measured close to the junction. Detailed 



analysis of this using the form of 6Gt given in Section 2.2|(i.e., neglecting the superfluid 



contribution here) leads to a mode occupation difference of [18], [19] 



2 



Evaluation of this quantity assuming momentum conservation during tunneling and 
using the parameters of the main tunnel junction in the Delft SQUID [7] leads to the 
following estimate for the tunneling contribution to the occupation number difference 
close to the junction 

33 < ANZ, < 43, (40) 

which is the same order of magnitude as the number of electrons changing modes in the 
bulk of the flux loop for this qubit (second row in Table [T|. Thus the tunneling dynamics 
of electrons close to the junction do not significantly alter the overall superposition size 
estimates made from microscopic consideration of the current carrying electrons in the 
bulk superconductor - both yield the same order of magnitude. Note, that while these 
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contributions are in general not simply additive, if the range of tunneling is restricted 
to the vicinity of the junction, however, then adding them does provides a reasonable 
bound on the total superposition size that consistently includes both bulk and tunneling 
electron contributions. 

4-3. Dirty Superconductors 

The analysis above has addressed ballistic superconductors at T = 0. We now 
consider the effect of impurities and defects on these calculations of cat size in terms 
of the number of participating microscopic degrees of freedom. First, we note that 
while superconductors may in general be subject to inelastic effects from magnetic 
impurities, their concentration in modern nanofabricated samples is, however, extremely 
low, to the extent that there are no impurity-induced states in the gap, the density 
of states is perfectly BCS-like [22] and the dephasing times of flux qubits [SU] are 
longer than expected in the presence of many magnetic impurities [31]. However, 
the superconductor may still be 'dirty' by virtue of elastic scattering of the Cooper 
pair electrons from impurities, which is likely due to the polycrystalline nature of an 
evaporated superconducting film. This results in typical diffusion constants of order 
D = 10~^ m^/s, corresponding to a mean free path i ~ 10~^ m, which may be shorter 
than the coherence length ~ hvp/A. We now consider how this scattering affects the 
mode occupation number difference (5nq in the dirty limit where i < (^o = 1-6 x 10^^ 
m in Al and 3.8 x 10"^ m in Nb). 

We need to analyze how the Green's function G(q, q') is modified by elastic 
scattering from impurities. We shall restrict our discussion here to the bulk contribution, 
GhvLik = Go+SGvg and not consider the tunneling contribution. Tunneling influences the 
Green's functions over a length corresponding to the appropriate dirty-limit coherence 
length ^o,D = a/^^o/S [32], which is necessarily smaller than in this limit. Thus, we 
expect a smaller tunneling contribution than in the ballistic case. For weak scattering, 
the modified bulk Green's function can be analyzed using the Dyson expansion 



where U is the electron-impurity interaction. Expanding as before in the small 
parameter |vs|/fcrit = IvsI^f^^/A leads to a zeroth order term GqUGq and a first 
order term GqUSGv^ + SGv^UGq. We take the electron-impurity interaction U to be 
defined by 



j 

where the sum goes over all impurities in the superconducting loop. The modified 
Green's functions may then be evaluated using the explicit solution for the zeroth order 



G — Gbuik + GhuikU Gbuik 



(41) 




(42) 
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Green's function, 



Go(k,r;k',r') = {2TTf5{k-k' 




{Ek(7, + Acr,)J < r - r' < +lh/3 

1 - ^ (Efecr, + Act,)] . -1^/3 < r - r' < 

(43) 

Carrying out the time and momentum space integrals, we find that the zeroth order term 
GqUGo is equal to zero and that the two contributions to the first order term cancel. 
Thus, to first order in the superfiuid velocity and the electron-impurity interaction, 
elastic scattering from impurities has no effect on the single electron Green's function 
and hence no effect on our estimate of the effective superposition size. 

While higher order terms in the Dyson expansion may also give rise to contributions 
that are first order in vs, a more complete analysis based on quasiclassical Green's 
function El] shows that these conclusions are independent of the impurity 

concentration as long as i ^ Xp, where Xp is the Fermi wavelength. The latter can 
always be assumed, because i c::^ Xp would put the material close to the Anderson 
metal-insulator transition in the normal state |35], hence the superconductor would 
not be a 'good metal' any more. The analysis proceeds as follows. In the dirty limit, 
the superfiuid velocity is related to the phase gradient by vs = -D (V0 + |f A), where 
D = ^vpi is the diffusion constant |3H|36]. The supercurrent density is still proportional 
to Vs and is found to be independent of D [31]. Furthermore, at a given current, the 
gauge-invariant phase gradient V0 + is proportional to 1/D. Now the angular 
dependence of the Green's function in the dirty limit is only weakly anisotropic, with 
a large s-wave component that does not depend on the direction of the circulating 
current (and hence does not contribute to the difference in mode occupation) and a 
small p-wave component that is proportional to i + [321 |33]. Using the above- 
mentioned dependence of the gauge invariant phase gradient on the diffusion coefficient, 
we then find that at a given current, the p-wave component of the Green's function is 
independent of the mean free path i. Futhermore, for the situation of a homogeneous 
circulating current, the angular dependence of the Green's function is restricted to s- 
and p-wave components only [33]. Thus the components of the Green's function relevant 
for evaluation of 6nq in superconducting fiux loops are independent of the impurity 
concentration and hence our estimates of effective superposition size apply also to the 
dirty superconductor regime. 



5. Discussion 



The results in Table[T]thus show that while not trivially small, the effective superposition 
sizes as bounded by the number of electrons that are in different modes in the two 
branches of the superposition are considerably smaller than previous estimates of the 
number of electrons carrying the supercurrent, which were based on simply counting 
all electrons within a London penetration depth of the surface [HI [71 |2l [T3|. In fact 
these estimates can now be replaced by the number of electrons in different modes 
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within the two branches, ANtot, since only this number of electrons is actually moving 
in opposite directions in the two branches. Specifically, our microscopic analysis of all 
electronic degrees of freedom shows that when electron indistinguishability is taken into 
account, the number of electrons responsible for the observed difference in supercurrent 
is actually far smaller than the number of electrons within the London penetration 
depth. Furthermore, there is a marked contrast between the fact that the values of 
AA^ are mesoscopic, while the two branches nevertheless have macroscopically distinct 
values of the two observables, persistent current and magnetic moment. We show below 
that this discrepancy is due to the Fermi statistics of the electrons. 

While the sizes obtained here are generally larger than typical estimates of size for 
superposition states realized in molecular and optical systems [U |6] and could reasonably 
be termed mesoscopic, they are still well short of anything that could be considered 
truly macroscopic. We therefore conclude that the superposition states in flux qubits 
are superpositions of relatively small numbers of particles that nevertheless result in a 
large difference in magnetic moment and current between the two branches because the 
electrons are circulating in opposite directions at high speeds (the Fermi velocity, which 
is of order lO^ms"^) around a path enclosing a relatively large surface area. The actual 
number of electrons that would be found to be behaving significantly differently in the 
two branches if one could measure them at the microscopic level is, however, seen to be 
quite modest. This situation is illustrated qualitatively in Fig. 1 by schematics of the 
shells at the Fermi surface that contain the occupied one-electron levels participating 
in Cooper pairing within each of the two branches of the superposition. The effect of 
realizing current fiows ±vs in the two branches is to shift each such Fermi shell by a 
distance ±mvg//i in k-space. Compared to mvp/h, the radius of the Fermi sphere, this 
distance is extremely small, so that the two Fermi shells overlap strongly. The vast 
majority of modes therefore have identical occupation numbers in both branches: only 
modes in the non-overlapping regions in Fig. 1 have a change in occupation number. 
However, these regions are on opposite sides of the Fermi sphere. Furthermore, simply 
stating that the two branches are separated by a slight shift of the entire Fermi sphere is 
incompatible with the fact that electrons are indistinguishable. Because of the electron 
indistinguishability and Fermi statistics, rather than shifting a large number of electrons 
by the small distance mw^/h, the shift is instead properly described as a relatively small 
number of electrons being moved all the way from one side of the sphere to the other 
and thereby changing their momentum by ~ 2/ci?. This is an extremely large velocity 
change because of the large diameter of the Fermi sphere, even though the small value 
of f s means that the volume of the non-overlapping regions and therefore the number 
of electrons moved are both small. This discrepancy between a large difference in the 
value of an observable quantity (current 5Ip or magnetic moment and the small 
number of particles actively involved in the superposition by distinguishably changing 
modes, ANtot, thus derives from the fermionic character of electrons. 

The corresponding gap between difference in observable quantities and number of 
particles involved might be smaller for large scale superposition states in a bosonic 
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I 

Figure 1. Schematic of tlie two sliells of single electron states at the Fermi surface 
involved in Cooper pairing for each of the two branches of the superposition jO) + |0). 
Each Fermi sphere is shifted by a distance ±2tovs /hin k-space and the shells are color 
coded as red and blue, with the overlapping region denoted by the color sum, purple. 
Note that the figure is not to scale, with the shift being exaggerated relative to the 
diameter of the Fermi spheres in order to be visible. Only modes in the non-overlapping 
regions (red and blue) have different occupation numbers in the two branches. 



system. Finding a physical system to realize this constitutes a challenge for future 
study, both experimental and theoretical. 

From Eq. (32) the maximum difference in occupation number for any one mode is 
bounded by /i |q ■ 5(vg) | /2A. This number is always small unless current differences 
are close to the critical current (more specifically, to the depairing current of the 
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Figure 2. A schematic of a large Josephson junction corrected to a narrow wire, 
maximizing the change in kinetic energy of the loop electrons and thus the ratio of 
AN/N. 

material), which is impossible in a system where the current passes through Josephson 
junctions. In fact, the ratio between critical and Josephson currents can be written as 
-^c,buik/-^c,j = 4:7ifRiy/Rs, where / is the superconducting condensate fraction (/ = 1 at 
T = 0), Rn the junction resistance in the normal state and Rs the Sharvin resistance of 
the superconducting material [22]. Given that Josephson junctions have high resistivity 
while superconducting metals have low resistivity, this ratio is usually large. One 
strategy to reduce it is to use large area Josephson junctions to obtain small Rn, 
connected by superconducting wires with small cross-section to maximize -R5, as shown 
in figure [2j Another option is to increase the number of modes available, by scaling 
up the physical dimensions of the system. However, the number of modes scales only 
linearly with system size. To reach truly macroscopic superposition states in this way 
would therefore require the physical dimensions of the flux qubits to be scaled up by 
many orders of magnitude. While this would certainly make it extremely challenging to 
maintain or observe any superposition behaviour in the presence of environmental noise, 
recent estimates indicate that coherence may still be observable in high inductance single 
Josephson junction loops with linear dimension of order 1 cm [37] . 

The results presented above show that the microscopic number of effective electron 
constituents in flux qubit superpositions with dimensions accessible to experiments 
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today, such that these constituents behave differently in the two superposition branches, 
is of order 10^ — 10^. This addresses a long-standing question related to macroscopic 
quantum coherence |2i |27j. We note that while this value derives from an operational 
measure that is independent of the physical system, asking only how many n- 
particle measurements are required to distinguish the two superposition branches, it is 
nevertheless useful to ask the further question of how these n-particle measurements may 
be realized and whether they may be decomposed into single particle measurements (for 
an example of the latter, see [IS])- Similarly, one may ask whether it is possible to control 
the AA^tot constituent degrees of freedom individually, as is assumed in most applications 
of a GHZ state, i.e., does the size of the superposition state reflect the number of useful 
dimensions (analogous to the distinction between useful and non-useful or "fluffy bunny" 
entanglement [38j)? Here we note that due to the large superconducting gap, actual 
realization of the distinguishability measure would require an n-electron Quantum Non- 
Demolition (QND) measurement that does not involve excitation across the gap. Such 
measurements have not been devised so far and would require bridging about one order 
of magnitude in energy between the qubit tunnel matrix elements and the energy gap. 
Nevertheless, the present bound on this measure does provide a microscopic analysis of 
the number of electronic constituents that could behave in a Schrodinger cat-like fashion, 
i.e., show (macroscopically) demonstrably different behavior in the two branches. It is 
then a further question as to whether and how this demonstrably different behavior of 
the individual constituents may be realized. 

6. Conclusion 

We have used the functional path integral formalism to connect the microscopic and 
macroscopic description of flux qubits and to derive expressions for correlation functions 
of creation and annihilation operators. This connection was then used to characterize 
the effective size of superposition states in these systems by evaluating the change in 
occupation number of electronic modes between the two branches of a flux superposition 
state that is characterized by macroscopic branch differences in physical observables. 
We showed that this quantity constitutes an upper bound on the distinguishability 
measure developed previously by two of us. The results obtained here for flux qubits of 
physical dimension accessible today show that the number of electrons, or equivalently, 
of Cooper pairs, that participate actively in the superposition behavior is of order 
10^ — 10^, considerably less than the total number of electrons participating in the 
supercurrent. This result shows that even if there is no intrinsic size or number scale 
limiting the existence of macroscopic quantum superpositions, the quantum statistics 
of the constituent particles can nevertheless be important for evaluating the effective 
number of particles participating in a superposition whose branches are characterized 
by macroscopically distinct observables. In particular, in the case of flux qubits, our 
analysis has revealed that it is a combination of the indistinguishability and Fermionic 
character of electrons, together with the large numerical value of the Fermi velocities. 
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that is responsible for the large change in current and magnetic moment per particle 
between branches, rather than a change in state of a macroscopic number of electrons. In 
summary, we see that a full microscopic treatment of the electrons in flux qubits shows 
that the superconducting flux qubit experiments performed to date provide neither 
a verification nor proof for the formation of true quantum superpositions on a scale 
beyond a few thousand microscopic constituent particles. The experimental quest for 
superpositions on the truly macroscopic scale as well as the verification or falsification 
of macrorealism by this route therefore remain open. 
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